Anomalous Spin Response in Non-centrosymmetric Compounds 
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We examine static spin susceptibilities x«/3(q) of spin components S a and Sp in the non- 
centrosymmetric tetragonal system. These show anomalous momentum dependences like Xxx(q) — 
Xyy(q) ~ Qx — an d Xxy(q) + Xyx(q) ~ (7a;<7y, which vanish in centrosymmetric systems. The 
magnitudes of the anomalous spin susceptibilities are enhanced by the on-site Coulomb interaction, 
especially, around an ordering wave vector. The significant and anomalous momentum dependences 
of these susceptibilities are explained by a group theoretical analysis. As the direct probe of the 
anomalous spin susceptibility, we propose a polarized neutron scattering experiment. 



Recently, the non-centrosymmetric heavy fermion su- 
perconductor attracts much attention. For a tetragonal 
compound CePtaSi, which has been intensively studied 
from both experimental and theoretical sides, the super- 
conducting transition has been observed at 0.75K be- 
low the antiferromagnetic transition temperature 2.2K at 
ambient pressure [1]. In addition, CeRhSi3 and CeIrSi3 
also show superconductivity by applying pressure 0, Q . 
Commonly, superconductivity in these compounds has 
been found around the antiferromagnetic phase like in 
centrosymmetric heavy electron systems []], 0, S Q • 

Theoretically, the non-centrosymmetric system is char- 
acterised by the Rashba-type effective spin-orbit cou- 
pling, which has an antisymmetric momentum depen- 
dence with respect to the spatial inversion [a, la]- As 
a characteristic feature of the model, it has been shown 
that in terms of a band splitting caused by the spin-orbit 
coupling, the uniform susceptibility has large Van-Vleck 
type contribution [7]. Furthermore, the thermal average 
of the spin operator of an electron with a momentum k 
does not necessarily vanish 0, 0] , although any magnetic 
moment disappears by canceling out spins in the para- 
magnetic state. The effect of electric spins on transport 
coefficients like spin Hall effect are intensively studied in 
the field of spintronics @, H, • 

For the normal state property of non-centrosymmetric 
heavy fermion superconductors, the difference from cen- 
trosymmetric systems by the Rashba-type spin-orbit cou- 
pling has not been suggested except for the quanti- 
ties mentioned above. Considering that in hydrostatic 
pressure, the superconductivity appears from the para- 
magnetic state with decreasing temperature, the spin 
fluctuation will be one of keys of superconductivity, 
and a normal state property characterising the non- 
centrosymmetric compound may relate with the mech- 
anism of the superconductivity. Furthermore, as well as 
transport coefficients, it is expected that the remaining 
electron spins affect the magnetic excitation in the non- 
centrosymmetric compound. 

In this Letter, we examine static spin susceptibil- 
ities Xa/3(q) of a- and /^-components of spin opera- 
tors, based on a simple Hubbard model including the 
Rashba-term. Significant momentum dependences like 



Xxx(<l) - Xyy(<l) ~ ql-Qy and Xxy(<l) + Xyxfa) ~ q x q y 
will be shown, where these susceptibilities vanish in the 
centrosymmetric system. In order to clarify the origin of 
momentum dependences, we carry out a group theoreti- 
cal analysis to show that the symmetry of momentum de- 
pendence of Xapio) i s identical with the representation of 
the product of spin operators included in a corresponding 
susceptibility. In order to observe the unusual momen- 
tum dependence of spin susceptibilities, we suggest a po- 
larized neutron scattering experiment, especially, around 
the magnetic instability. 

In the following, we describe a non-centrosymmetric 
system by the following Hamiltonian 



H — Hq + Hi , 



(i) 



H o = ^2 t( £k ~ + gk ■ ^W c L c k<r'» (2) 

\<L(T(T' 

#i =t/Y> iT n U) (3) 

i 

where c\^ a and c\. a are annihilation and creation opera- 
tors of an electron with a momentum k and a spin a. 
Here, and \i are the energy dispersion of electrons and 
the chemical potential, respectively, while gk describes 
the Rashba field satisfying g k =— g_ k , which breaks the 
inversion symmetry. Then, eigenenergies of the non- 
interacting part are given by £k± = £k ± |gk| — In 
addition, Hi is the on-site interaction term. 

In the non-interacting case, the electron Green's func- 
tion is given by 



G<°>(k,i 



gk 

l&l 



aGL 0) (k,i^ n ), (4) 



with G^ } (k,icj n ) = \ [(iu n - £ k +) 1 ± (i^n - £k-) 
where uo n is a fermionic Matsubara frequency. Here, 
since , \u n ) is expanded only in odd-power of |gk|, 

G^} (k, iuj n ) vanishes in the limit of |gk| — > 0. Summing 
up the Green's function with respect to the Matsub- 
ara frequency, the first term gives electron densities of 
(k, ±) -states, which determine shapes of two Fermi sur- 
faces, while the second term contributes electron spins 
of (k, ±)-states, whose directions are oposite from each 
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FIG. 1: Schematic views of charges (a) and spins (b) of elec- 
trons. Red and blue lines denote Fermi surfaces of £k- and 
Sk+, respectively. Red and blue wedges are spins of electrons 
belonging to £k- and £k+, respectively. Spins of £k- -electrons 
remain in the region between two Fermi surfaces. 



other due to the minus sign in front of {\uj n — £k-) 1 of 
G^(k, iu n ). These features of electrons are schemati- 
cally drawn in Figs. 1(a) and 1(b), respectively. From 
Fig. 1(b), it is understood that spins of electrons belong- 
ing to £k--band remain in the region between two Fermi 
surfaces, while a spin of an £k--band electron at a k can- 
cels out another spin of an £?k+-band electron at the same 
k within the inner Fermi surface. It is expected that the 
remaining electric spin affects the magnetic response. 

In order to examine the magnetic excitation in the non- 
centrosymmetic system, we consider the dynamical sus- 
ceptibility, defined as 



dr 



Afl n r 



<T T [(S?(r) - (0) - <S°»]>(5) 



where (• • •) means the thermal average of • • T r denotes 
the imaginary-time chronological ordering operator, and 
Q n is a bosonic Matsubara frequency. Here, charge and 
spin operators with a wave vector q are given by = 

\ Eka C L C k+qa and S% = \ £ kfl ^, G^, ^Cfc+q^ («=X, 

y, and z), respectively, where & a is an a-component 
of Pauli matrices. In centrosymmetric systems, all off- 
diagonal susceptibilities disappear. On the other hand, 
it should be noted that off-diagonal susceptibilities do 
not always vanish in non-centrosymmetric systems [ll), 
and even susceptibilities between spin and charge oper- 
ators Xac(q,i^n) and Xca(q,i^n) (a=x,y, and z) have 
non-zero values for Q n ^ 0. 

According to the diagramatic technique, we formulate 
X«/3(q, i^ n ) with use of the Green's function G^, (k, iw m ) 
given above. In the non-interacting case, the susceptibil- 
ity x^2(q, ifi n ) is calculated through a transformation 



from Xo 



t (q,in n ) defined by, 



Xcri<72Cr3<74 (q? l^n) 

~ T E E G °L ( k ' ^rn)G^ ai (k + q, iu m + ifi„){6) 



k m 



which corresponds to the one-bubble diagram [Toj]. In 
the interacting case with a finite U, the expression of 
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FIG. 2: (a) q-dependences of static susceptibilities, Xc?(q, 0), 

X&^q.O). xi°2(q,0), X^(q,0), and (xi°2(q,0)-xi° ) (q,0))/2, 
where insets show Fermi surfaces in the first quadrant of the 
Brillouin zone (left) and the path in q-space (right), (b) 
q-dependences of static susceptibilities, xi°J (q 5 0)? Xyz(^0), 
and x2c (q 5 0) with the momentum path shown in the inset. 



Xa/3(q, i^n) is given within the random phase approxi- 
mation (RPA) P, [H, EH as 



2x (0) (q,i^n)^l x (0) (q,i^n), (7) 



where matrix elements of U are given by U a p = 5a,pU aa 
with £7 CC = —U and [T^ = fT^ = [/^ = [/. 

Then, we calculate Xaa' (q? i^n) numerically, based 
on a two-dimensional non-centrosymmetric system 
with an energy dispersion £k = 2ti(cos/c ;E + 
cos ky) + 4t 2 cos k x cos k y and a Rashba- field gk = 
g(sin fc^, — sin fc^, 0). For parameters, £2/^1 =0.35 and 
g/ti=0.2 are chosen, we can reproduce quasi two dimen- 
tional Fermi surfaces of CePtaSi obtained by the band 
calculation [3, 15, [hJ. In Fig. 2, we show momen- 
tum dependences of static spin susceptibilities in the 
non-interacting case. Diagonal components of X^(q) 
are shown in Fig. 2(a), where Fermi surfaces and the 
momentum path are also shown in insets. Unlike in 
centrosymmetric systems, all momentum dependences 

of Xc?(q), xi?(q), XaS(q), and Xyy are different 
from each other. Especially, in a path from (0,7r) to 

(7r,0), Xxx (q) — Xyy (q) is antisymmetric around the mid- 
point (7r/2,7r/2), and it vanishes in a diagonal path from 
(7r,7r) to (0,0). Therefore, the momentum dependence 
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FIG. 3: Schematic views of particle-hole excitations. Red and 
blue lines denote Fermi surfaces of £k- and Sk+ , respectively. 
Red and blue wedges are spins of electrons belonging to £k- 
and £k+, respectively. In (b), particle- hole excitations with 
momenta q x and q y only contribute to Xyy fax) an d Xx°xfay), 
respectively, while every particle-hole excitaion is isotropic for 
the spin index a in (a). 



of Xxxfa) 



Xyy fa) nas the typical q£ — q'y symmetry. 

Likewise, momentum dependences of Xxyfa) + Xyxfa), 

(X { yz fa) - Xzy (q))/i, and (xi°i (q) - xi°i (q))/i are shown 
in Fig. 2(b), where the momentum path is depicted 

in the inset. In the figure, Xxyfa) + Xyx fa) is sym- 
metric around (0,0) and (27r,0), while it is antisymmet- 
ric around (7r,7r). Accordingly, the momentum depen- 



dence of Xxyfa) + Xyx } fa) is a q x q y -type. Similarly, it 
can be understood easily that momentum dependences of 

(X { yz (q) - Xzy fa))/' 1 and (x™ (q) - xi°i (q))/i are q y - and 
g^-types, respectively. Thus, spin susceptibilities have 
significant and unusual momentum dependeces related 
with corresponding spin indices. 

For an intuitive picture of anomalous momentum de- 
pendences of spin susceptibilities, we consider particle- 
hole excitations in the non-centrosymmetric system. Ac- 
cording to the prescription to calculate x^(q), contribu- 
tions to Xxx (q) and Xyy (q) are separated to two terms of 



,(°) 



convolutions, T (k, kc; m )G+ ; (k + q, iu m + iO n ) 

andTE m GL 0) (k,i^ m )GL 0) (k + q, \u m + iQ n ) , where the 
second part vanishes in the limit of the zero Rashba 
field. Here, we show schematic figures of these particle- 
hole excitations in Figs. 3(a) and 3(b), respectively. In 
Fig. 3(a), the mechanism of particle- hole excitations is 
not different from that of the centrosymmetric case. In 
this case, the contribution to Xalfa) is isotropic in the 
spin space. On the other hand, a particle-hole excita- 
tion with a momentum transfer q^ shown in Fig. 3(b) 
only contributes to Xyy fax) f° r low energy excitations: 
rotating the quantization axis to the — e y direction, the 
longitudinal spin excitation is only permitted, and other 
excitations are almost forbidden because of the energy 
loss from spin rotations against the Rashba-field. Simi- 
larly, if we consider a particle-hole excitation with a mo- 
mentum transfer q y shown in Fig. 3(b), only Xxx fay) 
has the contribution. This is consistent with a behavior 
in a path (0,0)-(0,7r) of Fig. 2(a), where the momen- 
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FIG. 4: q-dependences of static susceptibilities (a) Xzz{<±,0), 
(b) (x«-(q>0)+X»»(q,0))/2, and (c) (x^(q,0)-x TO (q,0))/2 
for several values of U. The position of the ordering wave 
vector is depicted in the inset of (b). 



turn dependence of Xxx fa) is considerably different from 
those of Xzzfa) and Xyy fa)- Considering these points, 
Xx°x fa) — Xyy fa) should have the q% — g^-type momentum 
dependence. Similar consideration will also work well for 

Xxyfa) +Xy°xfa). 

The above discussion is for the non-interacting case. 
By the interaction U magnitudes of anomaous Xxx fa) — 
Xyy fa) an d Xxyfa) +Xyxfa) should be changed. In order 
to investigate the effect of interaction on q-dependences 
of susceptibilities, we calculate Xapfa) within RPA using 
the same parameter set as in Fig. 2. With the critical in- 
teraction U c = 2.551ii, the paramagnetic state becomes 
unstable at an ordering wave vector Q = (0.8l7r, 0.257r). 
We show momentum dependences of Xzzfa), (Xxxfa) + 
X^(q))/2, and anomalous (Xxxfa) ~ Xyyfa))/^ in Figs. 
4(a), 4(b), and 4(c), respectively. Xzzfa) and (Xxxfa) + 
Xyy fa)) /2 are enhanced with increasing J7, but have dif- 
ferent q-dependences from each other. Furthermore, the 
amplitude of (Xxxfa) — Xyy fa))/ % increases considerably 
around the ordering wave vector Q in comparison with 
the corresponding value of the non-interacting case. In 
addition, we note that the amplitude of Xxyfa) + Xyxfa) 
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TABLE I: Classification of spin product S^S^ according to 
irreducible representations of the C4 V x AC group. Then, S c , 
S z , and {S x ,S y } belong to A+, A^~, and E~ irreducible rep- 
resentations, respectively. The third column shows the basis 
function in the momentum space. The superscript + (-) of ir- 
reducible representations expresses the even (odd) parity with 
respect to the time reversal. 



r 


S a S 13 


basis 




s c s c 


2 2 2 
Qz > Qx + Qy 




s z s z 






s x s x + s y s y 




A 2 + 


i(S c S z - S Z S C ) 


QxQy{Ql - qy) 


Bf 


s x s x - s y s y 


2 2 
Qx ~ Qy 


B 2 + 


s x s y + s y s x 


Q x q y 


E + 


{S y S z + S z S y , S Z S X + S X S Z } 


{q y qz,q z qx} 




{i(S c S x - S X S C ), i(S c S y - S y S c )} 




A 2 " 


s c s z + s z s c 


qxq y q z {ql - ql) 




i(S x S y - S y S x ) 




E" 


{i(S y S z - S z S y ), i(S z S x - S X S Z )} 


{qx,q y } 




{S C S* + S*S C , S c S y + S y S c } 





is also enhanced with increasing U. 

In order to clarify the origin of the unusual momen- 
tum dependence of the spin susceptibility, we carry out 
a group theoretical analysis. In the present case, the 
relevant group is x /C, where is the tetragonal 
point group without the inversion symmetry, and JC is 
the time-reversal symmetry group. From the definition, 
the momentum dependence of Xapi^ ifi n ) is determined 
by the expectation value of the right hand side. Noting 
that both spin and momentum are transformed by sym- 
metric operations of the group, the representation of the 
operator in the expectation value is T s <g> T m , where T s 
and T m are representations of the spin product and the 
momentum dependence, respectively. Since the operator 
giving a non-zero expectation value (O) belongs to the 
identity representation Ti (A+ in Table I) in the group, 
the relation T s (g) T m = Ti is obtained, so that finally 

r m = r s (8) 

is required. For the symmetry of the q-dependence of 
Xa/3(q, ifi n ), it is sufficient to know the representation 
of the spin-product S a S@, where the classification of the 
product is given in Table I. Therefore, representations of 
q-dependences of Xxx(<l) - Xyy(<l) and Xxy(<l) + Xyx(<l) 
are and Bj, respectively. These results are consistent 
with q-dependences of anomalous susceptibilities shown 
in Figs. 2 and 4. Thus, it is shown that the symmetry 
of the q-dependence of Xa/3(q, i^n) is identical with the 
representation of the spin- product S a S p . 

Finally, we suggest a polarized neutron scattering ex- 
periment to observe the novel anomalous momentum de- 



pendence of susceptibility in non-centrosymmetric heavy 
fermion systems, where even anomalous susceptibilities 
will be enhanced by the strong interaction as shown in 
Fig. 4. Among non-centrosymmetric heavy fermion com- 
pounds, the only available neutron scattering data are for 
magnetic structures of CePt3Si[17] and CeRhSis[18]. For 
CePtaSi, considering the ordering wave vector (0,0,7r), to 
observe the anomalous momentum dependence of sus- 
ceptibilities will be difficult, because the symmetry of 
q-dependence of every anomalous susceptibility is not 
the A^ representation. On the other hand, for CeRhSi3, 
whose ordering wave vector is (±0.437r,0,7r), it will be 
more promising to observe the anomalous momentum de- 
pendence of, especially, Xxx(q) — Xyy(°L) j us ^ above the 
magnetic transition temperature, since the susceptibil- 
ity will be enhanced around the ordering wave vector by 
approaching the magnetic transition temperature. 

In summary, we have examined the dynamical sus- 
ceptibility in the non-centrosymmetric tetragonal sys- 
tem. It has been shown that affected by the Rashba 
term, x^(q) - Xyy(^) and X^(q) + Xyxfa) have q 2 x - q 2 y - 
and g^g^-types of momentum dependences, respectively. 
The group theoretical analysis has been used to explain 
the peculiar feature in the non-centrosymmetric systems. 
Since the unusual susceptibility is also enhanced by the 
Coulomb interaction, it is desirable to observe the un- 
usual momentum dependence of susceptibility by the po- 
larized neutron scattering experiment. 

The author is grateful to P. Thalmeier for many valu- 
able discussions. 
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